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SUMMARY 

A study was made of the onset of flow oscillations occurring when a subcooled liquid 
undergoes a phase change under forced-flow conditions in a single-tube boiler. Experi- 
mental data were obtained with electrically heated tubes of variable length over a range 
of fluid velocities, temperatures, and pressures. The experiments were conducted at 
system pressures of 3 to 100 psia with degassed, demineralized water flowing vertically 
upward through the heated test tube. Visual studies were made in which boiling was 
simulated by injecting steam or air through the porous wall of a tube. A mathematical 
model predicting the necessary criterion for the onset of oscillations was formulated on 
the basis of local rates of bubble growth and collapse. Good agreement was obtained be- 
tween the experimental and theoretical results. 


INTRODUCTION 

The study of instabilities occurring during the change of phase of a flowing liquid has 
been the object of many investigations in recent years. Experimental data show that var- 
ious types of instabilities can result when boiling exists in a heated channel. High- 
frequency pressure fluctuations have been observed in subcooled boilers (ref. 1). These 
are believed to be caused by the vapor voids collapsing in the subcooled liquids. "Flow 
excursions,” first reported by Ledinegg (ref. 2), can result from improper matching of 
the pump and boiler hydraulic characteristics. An instability may also be reflected in a 
surface temperature excursion, which is no more than the well-known ’’burnout" phe- 
nomenon, generally defined as the transition from the nucleate to the film boiling regime. 
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The information presented herein was part of a thesis submitted by F. A. Jeglic to 
the University of Notre Dame in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy, October 1964. 



Another type of instability observed in boiling systems results in oscillatory flow behav- 
ior. Often the inception of flow oscillations determines the upper limit of the operating 
conditions of a given system. Consequently, much effort has been made to predict the 
onset of flow oscillations under various physical conditions. A review of the literature 
and state-of-the-art is given in references 3 and 4. 

The objective of this investigation was to study the onset of flow oscillations in a 
single-tube forced-flow subcooled boiler. For all cases considered, the fluid bulk tem- 
perature was below the saturation temperature at both the exit and the inlet of the boiler. 
Demineralized, degassed water was used as the boiling fluid. Test sections consisted 
of vertical electrically heated Inconel X tubes of various lengths. All experimental tests 
were conducted in the pressure range of 3. 0 to 100 psia with a liquid bulk velocity of 
2. 2 to 8. 3 feet per second for various inlet bulk subcoolings. Some special test sections 
were constructed which permitted flow observation and photography. For some of the 
visual studies, porous-wall tubing, which allowed boiling simulation by injecting steam 
or air through the porous wall, was used. Comparison between the results obtained with 
electrically heated test sections and the results of simulated boiling studies was made. 

A criterion predicting the necessary condition for the onset of flow oscillations in sub- 
cooled boiling is derived theoretically. The theory is based on the rate of growth and 
collapse of the vapor bubbles. Theoretical results were found to compare well with ex- 
perimental data of this investigation. 


APPARATUS 

General 

A schematic diagram of the closed-loop forced-flow-boiling heat-transfer facility is 
shown in figure 1 and is described in detail in reference 1. The system is pressurized 
or evacuated through the expansion tank by dry nitrogen gas or external vacuum pump, 
respectively. All component parts in contact with the test fluid were made of corrosion 
resistant materials to prevent fluid contamination and scale deposits. The fluid is circu- 
lated by a gear pump, which is powered by a variable speed drive. The power to the 
test section and the preheater is supplied by 250- and 70-kilowatt interchangeable 
alternating-current transformers, each regulated by a saturable reactor. The reactors 
provide fine continuous control of the power input. A counterflow heat exchanger serves 
as the system heat sink. No restriction at the inlet of the boiler was provided for the 
present series of runs. 
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Test Sections 


Electrically heated tubes . - The primary test section was constructed of 0. 020- inch- 
wall Inconel X tubing with a 0. 25-inch outside diameter. The tubing was especially 
rolled to limit the wall thickness deviation to less than 3 percent of the nominal wall 
thickness and had a maximum surface roughness of 20 microinches. The heated length of 
the test tube was 21 inches (100 diam) and was preceded by a 15-diameter unheated length 
to reduce the hydrodynamic entrance effects. Power leads were connected to copper 
flanges that were silver -soldered to the tube. High-temperature solder was used, and 
care was taken to ensure good contact and to avoid fillets in the heated portion of the test 
section. The flanges were bolted to the inlet and exit plenums and electrically insulated 
by Teflon and asbestos gaskets and seals. The exit plenum was fixed to the structural 
frame, while the inlet plenum was allowed to slide longitudinally to provide for thermal 
expansion. The test-section assembly was enclosed by a transparent shield to lessen 
external convective currents and provide protection in case of failure. 

A second test section of similar dimensions was built, but it was provided with a 
sliding copper flange. In this manner the length of the heated portion of the tube could be 

varied. 

Transparent test sections . - Three special 
test sections were built for visual studies. One 
was constructed of 0.21-inch-inside-diameter 
porous-wall stainless-steel tubing. A short glass 
tubing of the same inside diameter was fused to 
the metallic tube at the exit. The porous -wall 
tube, with 5- micron-diameter pores, was en- 
closed by a 1-inch-diameter tube, which formed 
a chamber connected to a degassed steam or air 
supply. This test section is shown in figure 2. 

Two annular test sections were constructed; 
the inside tubing was metallic and the outside 
tubing was glass. In the one shown in figure 3, 
the inside tube was 0. 010-inch- wall, 1/4-inch- 
outside- diameter Inconel X tubing connected to an 
electric power supply; in the other, shown in 
figure 4, the inside tube was 1/4-inch-outside- 
diameter stainless-steel tubing with a porous 
wall (5-micron-diam. pores) through which steam 
could be injected. For both test sections, the 

Figure 2. - Tubular, porous-wail test section with trans- inside diameter of the glass tube was 5/8 inch, 

parent extension (used in photographic studies). 



From inlet plenum 


*0-71529 
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Figure 3. - Electrically heated, transparent annular test Figure 4. - Annular, transparent test section with porous 

section (used in photographic studies). inside wall (used in photographic studies). 


Instrumentation 

The average values of pressures at various locations in the facility were measured 
with a number of Bourdon tube gages. Fluid bulk temperatures were measured by 
copper-constantan thermocouples and recorded on a multichannel self-balancing potenti- 
ometer. A number of 36-gage Chromel-Alumel thermocouples, spotwelded to the outside 
surface of the boiler tube, were used to record the surface temperature. A thermocouple 
spotwelded 1/ 4 inch from the exit end of the heated length was connected to a separate 
recording instrument and served as an overtemperature control; that is, an unexpected 
large temperature excursion would automatically deenergize the power supplies. A 
dynamometer-type wattmeter was selected to measure the power input to the test section. 
A multirange, electronic tube voltmeter was used to measure the voltage drop across 
the tube. The average flow rate was determined by the signal generated by a turbine- 
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type flowmeter and read on the flow converter and a digital frequency counter. A de- 
tailed description of the instrumentation is given in reference 4. 

Strain- gage-type transducers with an accuracy of 1/4 percent of full scale were used 
to measure the instantaneous pressures at various locations of the loop (see fig. 1, p. 3). 
The signals were continuously recorded on a multichannel oscillograph. The flow tran- 
sients were obtained by connecting the signal, generated by the turbine flow meter, to 
the oscillograph. The response time of the flowmeter was 3 to 5 milliseconds, which is 
sufficient to obtain the instantaneous values of flow over the range of conditions studied 
herein. One of the thermocouples spotwelded to the wall was also connected to the oscil- 
lograph and provided qualitative measurement of the surface temperature variation. 
Whereas the response of the thermocouple was sufficient to detect the temperature oscil- 
lations, the recorded signal was somewhat obscured by the alternating-current pickup 
from the power supplies. 

A high-speed camera, capable of speeds up to 6000 frames per second, was used to 
photograph the transparent test section during the run in the visual studies. A photocell 
was mounted next to the test section and its signal recorded on the oscillograph. A light 
source could thus be recorded simultaneously on the film and the oscillograph. This 
enabled synchronization of the high-speed movies with the flow and pressure traces. 


PROCEDURE 

Preparations 

Prior to filling, the system was flushed with deionized water, purged with nitrogen, 
and evacuated. It was then filled with deionized, deaerated water. The liquid was cir- 
culated through the system and further degassed by boiling in the test section and venting 
to the atmosphere. Water samples were taken and analyzed on a gas analyzer for gas 
content. Gas content of 3 ppm (by weight) was considered a maximum permissible 
amount; the degassing process was continued until the air content was well below this 
value. 

The electric resistance of the test-section tubing was measured a priori, and the 
values obtained agreed well with the electric resistivity given in reference 5. It was also 
found that for the range of temperatures of this investigation the resistance may be con- 
sidered independent of temperature. 
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Calibrations 


All instruments were calibrated prior to installation and were periodically checked. 
The temperature recording instruments were calibrated with a standard potentiometer. 

All pressure gages were checked on the deadweight tester. The transducers were me- 
chanically and electrically balanced and checked against a known pressure. Under 
steady-state conditions, a good correspondence was obtained between the readings of the 
pressure gages and those recorded on the oscillograph. The flowmeters were calibrated 
prior to their installation. The flow converter and the digital frequency counter were 
calibrated by a standard signal generator. The wattmeter and the voltmeter, both used 
for input power measurement, were calibrated with a known source. Their rms charac- 
teristics were checked with known square, triangular, and sinusoidal signals. 

The power input calculated from the wattmeter measurements and from the voltmeter 
reading (together with the known electric resistance of the tube) agreed well in all cases 
in which both readings were taken. These values were also in good agreement with the 
calculated enthalpy increase of the fluid, which served as a further check on the power, 
flow, and bulk temperature readings. When the liquid was preheated and all the heat 
load removed from the test section, it was possible to check the wall temperature and 
fluid bulk temperature instrumentation. In effect, therefore, all pertinent instrumenta- 
tion was continuously checked. 


Experimental Technique 

E lectrically heated tubes . - The normal operating procedure consisted of setting the 
flow rate, inlet bulk temperature, and exit pressure to the desired values and then slowly 
increasing heat flux to the boiler until the oscillations set in or net vapor was generated, 
whichever occurred sooner. Flow rate and exit pressure were adjusted during the run 
to maintain the preset value. For the runs with constant amplitude oscillations the sys- 
tem was permitted to oscillate for a substantial period of time before the heat flux was 
reduced. The runs with increasing amplitude oscillations, however, were terminated by 
the automatic power shutoff, which was energized by the surface temperature excursion. 
Any one of the three controlled variables, flow rate, inlet bulk temperature, or exit 
pressure, was then changed and a new run begun. All recording instruments were turned 
on before the onset of oscillations and left on for the duration of the run. 

Once the heat flux corresponding to the incipience of oscillations was determined for 
a given set of conditions, the mode of operation was changed for a few runs to determine 
its effect on experimental results. Either the flow rate or pressure was raised initially, 
and the heat flux was set at the value previously found to induce oscillations, whereupon 
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TABLE I. - EXPERIMENTAL DATA WITH CONSTANT LENGTH, 
ELECTRICALLY HEATED TUBE 
[Tube inside diameter, 0. 21 in. ; length, 21. 0 in.] 


Run 

Inlet bulk 

Exit bulk 

Exit pressure, 

Inlet liquid 

Heat flux (at onset 


temperature, 

temperature, 

psia 

velocity, 

of oscillations), 


°F 

°F 

(a) 

ft/sec 

Btu/(sec)(sq ft) 

b 20 - 1-1 

79 

222 

20.0 

2.20 

49. 1 

b 20 - 1-2 

148 

223 



2.20 

34. 1 

c 20 - 1-3 

175 

— 



2.20 

— 

b 20 - 2-1 

68 

211 



4. 00 

93. 5 

b 20 - 2- la 

67 

208 



4. 00 

91.7 

b 20 - 2-2 

94 

216 



4. 00 

79. 0 

b 20 - 2-3 

115 

219 



4.00 

68. 5 

b 20 - 2-4 

146 

— 



4.00 

55.8 

c 20 - 2-5 

183 

— 



4. 00 

— 

d 20 - 3-1 

61 

211 



6.25 

146.0 

d 20 - 3-2 

91 

220 



6.25 

129 

d 20 - 3-3 

108 

222 



6.25 

113 

d 20 - 3-4 

130 

— 



6.25 

94.7 

e 20 - 3-5 

150 

— 



6.25 

87.6 

c 20 - 3-6 

182 

--- 



6.25 


d 20 - 4-1 

64 

219 



8.25 

201 

d 20 - 4-2 

99 

217 



8.25 

161 

f 20 - 4-3 

124 

— 


( 

8.25 


g 40 - 1-1 

80 

— 

40.0 

2.20 

75.6 

c 40 - 1-2 

152 

— 

40.0 

2.20 


b 40 - 2-1 

70 

266 

40.0 

4. 00 

125 

c 40 - 2-2 

100 

— 

40.0 

4.00 


*V - 1-1 

80 

140 

3.0 

2.20 

19.7 

V - 2-1 

71 

140 

3.0 

4.00 

39.0 

d V - 3-1 

70 

140 

3. 0 

6.25 

62. 5 

d V - 4-1 

70 

140 

3.0 

8.25 

85. 6 


a For all runs at 60 and 100 psia exit pressure, over the same range of liquid veloc- 
ities and maximum inlet subcooling, net quality was easily obtained without ob- 
serving any oscillatory behavior in the subcooled regime. 

D Constant- amplitude oscillations inflow, pressures, and surface temperature. 
c Net quality vapor easily obtained. 

^Oscillations with increasing amplitude that resulted in temperature excursion. 
e Irregular oscillations apparent; net quality easily obtained. 

f 

No oscillations detected in subcooled regime; net quality easily obtained. 
^Irregular oscillations apparent; net quality easily obtained. 



the flow or pressure was reduced slowly until the system was rendered unstable. It was 
found that the onset of oscillations was independent of the mode of operation for the range 
of this investigation; consequently, the instability was obtained by increasing the heat 
flux for all subsequent testing. 

A similar procedure was followed for the test section of variable length, except that 
the inlet temperature and exit pressure were not varied. Instead, the length of the test 
section was changed between the runs. 

Visual studies . - As before, the flow rate, inlet bulk temperature, and exit pressure 
were set to the desired values. With electrically heated test sections, the heat input was 
slowly increased until oscillations were recorded on the oscillograph. The high-speed 
movie camera was then activated as well as the light source. Most of the runs were 
taken at 1000 frames per second. For some runs the heat flux was reduced during the 
run in order to observe the stabilization of the system. 

The same variables were controlled with test sections in which two-phase flow was 
obtained by steam injection through the porous tube wall. In these studies, steam or air 
flow rate was increased slowly instead of electric power. The steam flow rate was either 
kept constant, decreased, or increased during the run, depending on whether it was de- 
sired to observe the neutral oscillations, the development of the oscillations, or the de- 
struction of the oscillations, respectively. 


EXPERIMENTAL RESULTS 
Electrically Heated Tubes 

The experimental results obtained with a constant length, electrically heated tube 
are given in table I. Each series of runs consisted of a number of runs at constant exit 
pressure and inlet bulk velocity, but at decreasing inlet subcooling. When the inlet sub- 
cooling was such as not to cause the flow to oscillate before net quality was being gener- 
ated, the series was terminated; the preheater load was then reduced to zero, the exit 
pressure or fluid mass rate was fixed at a new value, and a new series was begun. The 
values given in table I correspond to the conditions at the onset of oscillations, when 
these occurred in the subcooled boiling regime. 

As shown in table I, apparently different types of oscillations were recorded in sub- 
cooled boiling: 

(1) During operation at steady-state conditions, an incremental increase in heat flux 
induced flow, pressure, and surface temperature oscillations. These developed quickly 
and reached maximum amplitude in a few cycles. The oscillatory condition could then be 
maintained indefinitely. Typical traces are shown in figure 5, where the variations of 
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(a) Run 20-1-2. (b) Run V-l-1. 

Figure 5. - Typical traces illustrating constant-amplitude oscillations. 

pressures and flow are given as functions of time. Further increase in heat flux would 
result in either surface temperature excursion or net vapor generation. 

(2) During operation at steady-state conditions, an incremental change in heat flux 
induced oscillations with increasing amplitude and would result in surface temperature 
excursion. An example of the surface temperature, pressures, and flow variation with 
time is shown in figure 6. It is believed that these oscillations would reach a constant 
amplitude if the heat capacitance of the wall could be increased so that the excursion 
could be decreased sufficiently to prevent the power trip. 

(3) Irregular flow oscillations occurred, but their inception is difficult to define be- 
cause of their relatively small amplitude. A further increase in heat flux would result in 
net vapor generation. 

The experimental results show that for the range of this investigation the onset of 
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Figure 6. - Typical traces illustrating oscillations with increasing amplitude and resulting in surface temperature excursion. 


flow oscillations in forced-flow subcooled boiling is a strong function of inlet subcooling 
and system pressure. An increase in inlet subcooling or a decrease in system pressure 
tends to induce subcooled boiling oscillations. It was found that when the inlet bulk tem- 
perature and/or exit pressure reached some critical value the oscillations could not be 
obtained in subcooled boiling regime; under such conditions net quality vapor could be 
generated at the exit. It is quite possible that in boiling studies under elevated pressures 
this minimum required pressure was exceeded, and consequently the subcooled boiling 
oscillations were not observed. 

The data obtained with a variable length, electrically heated tube are given in 
table n. For these runs the exit pressure and inlet temperature were held constant, and 
the length of the test section was varied. The results indicate that the decrease in the 
boiler tube length tends to decrease the exit bulk temperature at which the flow oscilla- 
tions will develop. It should be noted that, whereas the total enthalpy rise of the fluid 
decreases with decrease in length, the heat flux actually increases. 
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TABLE n. - EXPERIMENTAL DATA WITH VARIABLE 


LENGTH ELECTRICALLY HEATED TUBE 


[Tube inside diameter, 0. 21 in. ; exit pressure, 20. 0 psia; 
inlet liquid velocity, 2. 20 ft/sec.] 


Run 

Test-section 

length, 

in. 

Inlet bulk 
temper- 
ature, 

°F 

Exit bulk 
temper- 
ature, 
°F 

Heat flux (at onset 
of oscillations), 
Btu/(sec)(sq ft) 

1-1 

18.25 

72 

216 

60.8 

1-2 

14.0 

70 

211 

76. 8 

1-3 

10. 5 

72 

214 

98.2 

1-4 

10.0 

70 

208 

104 

1-5 

7.0 

70 

199 

147 

1-6 

4.0 

72 

189 

234 


Visual Studies 

Initially, for visual studies, the porous -wall tube with a glass extension was used. 
Liquid flow rate, pressure, and inlet temperature were adjusted to values for which data 
had been obtained with the electrically heated test section. The rate of steam injection, 
through the porous walls of the tube, was then slowly increased until flow oscillations 
were recorded on the oscillograph, whereupon the movie camera was activated. The 
steam injection rate was not measured, but it was possible to calculate the heat input by 
the enthalpy increase of the fluid. The correspondence of conditions with the two test 
sections was excellent; not only did the oscillations occur at the same values, but also 
the traces were similar. A comparison of the recorded traces is illustrated in figure 7. 
The flow regime was observed to be bubbly at steady-state conditions, just preceeding 
the onset of oscillations. As a small increase in steam injection rate resulted in record- 
ing oscillations on the oscillograph, a vapor slug appeared in the transparent portion of 
the test section. The liquid-vapor interface along the wall at the beginning of the slug 
was quite wavy and irregular. As the slug proceeded up the tube, the liquid film on the 
wall became increasingly thinner with the interface continually growing smoother. The 
vapor slug ended abruptly and was followed by bubbly two-phase flow. Immediately be- 
hind the vapor slug the bubbles were small and few but slowly grew in size and number 
until a new slug appeared. The length and the frequency of the slugs depended on the ex- 
perimental conditions. However, for all cases the frequency of vapor slug formation 
corresponded with the frequency of the pressure and flow oscillations recorded on the 
oscillograph. A sample trace is shown in figure 8 (p. 14). The effect of steam- 
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Flow increase 

Temperature and pressure increase 

(a) Steam injection. (b) Electrical heating. 

Figure 7. - Comparison of traces obtained with electrically heated test section and with steam injection. 
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■+■ Flow increase 

Temperature and pressure increase 

Figure 8. - Typical trace obtained with tubular, porous wall test section. 
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injection-rate decrease and subse- 
quent increase is illustrated in fig- 
ure 9. It can be seen that an in- 
crease in steam injection rate can 
cause a rapid development of oscilla- 
tions, while decreasing the steam 
rate results in a stable system. 

The high-speed movies were 
carefully examined and compared to 
the corresponding traces. From 
such a comparison it was concluded 
that, for every cycle, the bubbly 
flow regime existed when the flow 
was at the maximum. The flow began 
to decrease just before the head of 
the vapor slug appeared in the trans- 
parent portion of the tube, continued 
to decrease as the slug grew and 
progressed toward the exit plenum, 
and reached a minimum before the 
appearance of the tail end of the 
vapor slug. A complete cycle is 
shown in figure 10, where the pic- 
tures, obtained with the high-speed 
movie camera, are superimposed on 
the trace. It should be noted that 
figure 10 represents one cycle of the 
run given in figure 8. The photo- 
graphs represent approximately a 
1-inch length of the transparent por- 
tion of the tube, which is immediately downstream of the simulated boiling portion of the 
tube. Even though these movies do not show the actual development of the vapor slugs, 
they verify that flow oscillations and the ” slugging” appeared simultaneously. 

In a few runs, air instead of steam was injected through the porous wall of the tube. 
This represented an adiabatic and nearly an isothermal two-phase flow. The most favor- 
able conditions for oscillatory flow were selected, yet over the full range of air injection 
rate the system remained stable. As the air injection rate was increased, the flow 
regime changed from bubbly to wavy annular to smooth annular. No slugging flow pattern 
could be obtained. This was an important observation which led to the conclusion that a 




P in * 20 psi/in. 


1 in. 
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rF = ~45 {lb m /hr)/in. 


— Flow increase 

Temperature and pressure increase 


Figure 9. - Effect of steam-injection-rate increase and decrease on 
stability in simulated subcooled boiling. 
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two-component two-phase system will not result in oscillatory behavior in the range of 
this investigation. Thus, the analogy between single-component and two-component two- 
phase systems, so often employed, could lead to erroneous conclusions. 

The fluid flow pattern behavior at the inception of oscillations in subcooled boiling 
was visually examined in the two annular test sections, one of which was electrically 
heated while the other used steam injection through a porous wall to simulate boiling (see 
figs. 3 and 4 (p. 5)). The preliminary results obtained by the two methods agreed quali- 
tatively as well as quantitatively. It was observed visually that, at the inception of oscil- 
lations, there is an abrupt transition from bubbly to annular flow at the exit of the test 
section. This is illustrated in figure 11. The transition is so sudden that an interface 
between the two flow patterns can easily be defined, the annular flow portion giving the 
appearance of a vapor slug. This '’interface’' oscillates up and down; the frequency of 
oscillation agrees well with the frequency recorded on the oscillograph. 


ANALYSIS 

When a subcooled liquid is brought in contact with a heating surface, certain temper- 
ature gradients will be established within it. If the heating surface temperature is above 
the saturation temperature, vapor voids may be formed at the surface. These voids may 
be ’’torn away” from the surface by one or a combination of the existing forces and con- 
dense in the cold liquid. 

When a subcooled boiling liquid flowing through a heated channel is considered, the 
radial temperature profile through the fluid will in part determine how the vapor voids 
formed at the wall will grow or diminish in size as they penetrate into the lower tempera- 
ture fluid away from the heated surface. If the conditions are such that the vapor voids 
emitted from the wall quickly decrease in size and finally collapse, the two-phase flow 
pattern will give the appearance of bubbly flow. However, if the vapor voids cannot col- 
lapse in the cold liquid for one reason or another, it is feasible for these voids (or bub- 
bles) to agglomerate and form larger voids, until they finally give the appearance of a 
vapor slug at some point in the boiling channel. 

As discussed previously, it was observed experimentally that in subcooled boiling 
the onset of oscillations was simultaneous with the appearance of the slug flow pattern in 
the boiler tube. When the flow was oscillatory, an oscillating well-defined "interface” 
between different flow regimes was observed. 

Thus, it is postulated that for oscillatory flow behavior it is necessary for the radial 
temperature distribution in some section of the boiler tube to be such as to permit the 
bubbles formed at the heated surface to grow. Obviously, the bubbles in the subcooled 
liquid at the center of the tube will tend to decrease in size, while those in the high- 
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temperature boundary layer will tend to grow. If the combined effect of the boundary 
layer and the cold liquid outside the boundary layer results in net growth of the vapor 
bubbles, then and only then is it possible for the bubbles to agglomerate and form a vapor 
slug. The formation of the slug would in turn satisfy the necessary condition for the in- 
ception of oscillations. Therefore, if the necessary criterion for the inception of oscil- 
lations in subcooled boiling is to be found, it becomes imperative to determine the radial 
temperature distribution and examine its effect on the bubble growth. 

At the inception of subcooled boiling, the vapor voids formed at the heated wall are 
small and few because of the relatively thin thermal boundary layer and the high liquid 
bulk subcooling. The bubbles that escape from the heated wall quickly collapse in the 
cold liquid. It would then be expected that in this region the bubbles would have only a 
small effect on the hydrodynamic behavior of the fluid. Experimental evidence shows 
that the bubbles remain small for some time after the inception of subcooled boiling, even 
at very low pressures. As the bulk temperature of the fluid is further increased down- 
stream in the tube, a point is reached where the bubble growth becomes of great impor- 
tance. This region is developed quickly, as can also be deduced from the theoretical 
bubble growth analysis, which proves the bubble growth to be an exponential function. 
Whereas the hydrodynamic and bubble growth phenomena are in reality coupled, it may 
be postulated that initially the liquid hydrodynamic phenomena are of prime importance, 
while the bubble growth effects control the fluid behavior in the last stages of the sub- 
cooled boiling regime. 

Thus, when the fluid behavior is strongly dependent on the bubble growth across the 
tube, the combined effect of the boundary layer and the cold liquid outside the boundary 
layer results in some net growth of the vapor bubbles. This condition, however, was 
postulated to be the necessary condition for vapor slug formation, that is, for inception 
of oscillations. Then, if the transition from the hydrodynamic to the bubble growth re- 
gion is relatively short, the thermal behavior of the hydrodynamic region in which the 
bubble growth has a negligible effect establishes the necessary condition for inception of 
oscillations. Therefore, in the analysis to follow, the temperature distribution across 
the tube was determined by neglecting the effect of bubbles; that is, the fluid was as- 
sumed to be in the liquid phase. After the temperature distribution was obtained, its 
effect on bubble growth was examined. 


Temperature Profile 


The problem of heat transfer across turbulent incompressible boundary layers has 
been dealt with by many authors. A discussion of various analyses and their limita- 
tions is presented in reference 6. An analysis for obtaining velocity and temperature 
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distributions in a tube with turbulent flow that applies well for water and takes into ac- 
count the variation of viscosity with temperature is the investigation of Deissler (refs. 7 
to 10). The results of Deissler, which are applicable only for fully developed turbulent 
flow, were used in the present work to predict the temperature profiles. Deissler's 
analysis was used primarily for convenience and because of its applicability to water and 
the present experimental conditions. Its development is briefly reviewed in appendix B. 

As shown in appendix B, the temperature profile can be calculated from the equa- 
tion (B7), which can be rewritten as 


t = (1 - 


(1) 


where 


t t = 41. 0 


^ (N Re r°- 8 (N E x 0 ' 4 


Pr' 


+ tv 


( 2 ) 


(All symbols are defined in appendix A. ) Equation (2) is the well known Dittus-Boelter 
equation, which is applicable to fully developed turbulent flow through tubes. 


Bubble Growth 

Determining the growth of the vapor bubble has been the objective of many investiga- 
tions in recent years. Scriven (ref. 11), Forster and Zuber (ref. 12), and Plesset and 
Zwick (ref. 13) have all obtained solutions to the problem, all of which are in good 
agreement, even though the method of analysis was different. Very recently, Yang and 
Clark (ref. 14) arrived at the following simple solution for bubble growth constant £ by 
the source theory: 


f = (t - t ) (3 

2 'V'fg 

where Pp p y , c , h^, and t g are densities of liquid and vapor, specific heat, latent 
heat of evaporation, and interface temperature, respectively. Equation (3) was derived 
for the asymptotic stage of the bubble life, which is the time of interest in the present 
analysis. The results of Yang and Clark are in good agreement with other theories, as 
discussed in reference 14. As shown in references 11 to 13, the bubble growth in the 
asymptotic stage may be expressed as 
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R(0) = 


(4) 


where a ^ is thermal diffusivity and 9 is time. Differentiation of equation (4) yields 
the rate of change in the bubble radius: 


R 


-V? 


Substituting equation (3) in equation (5) gives 


R = ^I (t - t ) 

2 p..h» 10 s 


Vfg 


The average bubble growth rate over any cross section of the tube is then 


(5) 


(6) 


R 


av 



(7) 


Substituting equation (6) in equation (7) and assuming the properties to be constant give 


R av 


_ vgr pfy, i 

r?>/e p v h fg 



( 8 ) 


Equation (8) is an expression for the average rate of bubble growth (or collapse) over a 

cross-sectional area of the tube in which temperature gradients exist. It is evident that 

if the liquid temperature is greater than the interface temperature t g the bubbles 

will grow; if t < t , the bubbles will tend to collapse, 
s 


Necessary Stability Criterion 

As discussed previously, the visual experimental studies showed that, in subcooled 
boiling over the range of this investigation, the flow oscillations resulted simultaneously 
with the appearance of the "slugging" flow pattern in the boiling tube. This observation 
led to the assumption that the slugging was the necessary condition for onset of flow 
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oscillations. It was further postulated that slug formation is only feasible when the con- 
ditions at some cross section of the tube are such as to yield some net growth of the vapor 
voids formed at the heated surface. This is equivalent to saying that the right side of 
equation (8) is greater than 0, or in the limit is equal to 0. Mathematically this may be 
expressed as 


f Vi (t - t )r dr = 0 (9) 

• / 0 

where t is the fluid temperature defined by equation (1) and t„ is the bubble interface 

ij 

temperature. Equation (9) is then the necessary condition for the onset of flow oscilla- 
tions for subcooled liquid boiling in a round tube. 


Calculation of Necessary Stability Criterion 

The procedure in determining the conditions under which the necessary criterion for 
stability (eq. (9)) would be met was as follows. 

The physical dimensions of the tube, the fluid inlet velocity, and the system pressure 
were selected such as to correspond to the experimental data. Inlet bulk temperature 
and heat flux were then arbitrarily chosen. From these values the bulk temperature at 
the exit was calculated from the heat balance: 


^b, ex in + 


2q i L 

: / r i u b 


( 10 ) 


For uniform heat flux it is obvious that the most critical case in the tube is at the exit, 
where the fluid is at the highest bulk temperature. If the necessary stability criterion is 
to be satisfied at any cross section of the tube, it will first be met at the exit. There- 
fore, exit conditions were of prime consideration. 

Based on the exit bulk temperature, the surface temperature was calculated from 
equation (2). The initial value of shear stress at the wall was estimated from 1/7 power 
law, which implemented the calculation of the initial value of the heat-transfer parameter 
from 


0 = 



C T.t. 
pi 1 


( 11 ) 
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For y* < 26, the velocity and temperature parameters were computed as a function of 
y* by simultaneously iterating equations (B3) and (B4). As suggested in reference 10, a 
value of n = 0. 124 was used. For water, for the range of interest, it was found the con- 
stant m in equations (B3) and (B4) was m = -1. For 26 < y* < r*, equations (B5) 
and (B6) were used to evaluate the velocity and temperature parameters, when K = 0. 36, 
as recommended in reference 10. 

The universal temperature and velocity distributions were thus obtained over the full 
range of y* for the estimated value of t^. These calculated values were then used in 
equation (B9), the result of which made it possible to determine a new value of Tj by 
equation (B8). This new value of shear stress at the wall was used to recalculate the 
heat-transfer parameter, and the process was repeated until the assumed and calculated 
shear-stress values were nearly equal. The temperature parameter values were then 
converted into temperatures by equation (1) and matched with corresponding values of r 
by 


r 


r i 


Mi 




( 12 ) 


The temperature as a function of radius was thus found for a given set of conditions. Be- 
fore the integral of equation (9) could be evaluated, however, the interface temperature 
t had to be known. This value was assumed to be the local saturation temperature. 
Since the pressure was considered constant across the radius in the derivation of equa- 
tion (9), the saturation temperature was also considered independent of the radius. 

Hence, with the calculated temperature profile and the assigned interface temperature it 
was possible to evaluate the integral in equation (9). 

A new value of heat flux q^ was chosen, and the cycle was repeated until the value 
of the integral changed the sign. The values of the integral thus obtained were plotted 
against the corresponding values of heat flux, and the result was a value of heat flux at 
which the integral is equal to 0. 

All computations were performed on an IBM 7094 digital computer. Simpson's rule 
was used for all integrations, with a minimum number of increments being regulated by 
a maximum allowable error of 0. 1 percent. All convergence tests were based on 0. 01 
percent error. The program was debugged and checked by results obtained independently 
on an IBM 610 computer. The range of variables used in the calculations was limited to 
the range of experimental data. The variables used included pressure, velocity, bulk 
temperature, and boiler length. The results are discussed in the next section. 
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Heat flux, qj, Btu/(sec)(sq ft) 


DISCUSSION OF RESULTS 


Effect of Inlet Bulk Temperature 



Inlet bulk temperature, °F 


(b) Bulk velocity, 4.00 feet per second. 

Figure 12. - Comparison of experimentally determined and theoreti- 
cally predicted necessary heat flux for onset of oscillations as 
function of inlet bulk temperature. Exit plenum pressure, 20 psia. 


In calculating the heat flux that 
satisfies the necessary criterion for 
stability (eq. (9)), inlet bulk tempera- 
ture was taken first as the independent 
variable. The computed necessary 
heat flux is plotted against inlet bulk 
temperatures in figure 12 for four 
different velocities. The experimen- 
tal data, given in table I (p. 8), are 
superimposed on the graphs for com- 
parison purposes. Good agreement 
between experimental and analytical 
results is obtained, except at high 
inlet bulk temperatures. This means 
that over the full length of the tube the 
bulk temperature of the fluid in this 
case is relatively closer to the satu- 
ration temperature. Consequently, 
the conditions are more favorable for 
formation of larger and more numer- 
ous vapor voids at the heated surface. 
Thus, in this region the bubbles be- 
come of increasing importance, which 
violates the single-phase assumption 
postulated in the analysis. 

Furthermore, when the inlet sub- 
cooling is small, the heat input nec- 
essary for net vapor generation at the 
exit is much lower than that for large 
subcooling. Therefore, the axial rate 
of change of void fraction would be 
expected to be less, resulting in a 
transition from bubbly to annular flow 
pattern, as depicted in figure 13(a). 

It is speculated, however, that if the 
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Heat flux, q f , Btu/(sec)fsq ft) 


p 



60 80 100 120 140 160 180 

Inlet bulk temperature, °F 

(d) Bulk velocity, 8.25 feet per second. 


Figure 12. - Concluded. 



(a) Unfavorable. (b) Favorable. 

Figure 13. - Flow pattern development for inception 
of oscillations. 


flow oscillations are to occur in sub- 
cooled boiling, this transition must be 
abrupt, as illustrated in figure 13(b). 
This means that if the flow oscillations 
are to occur, the axial rate of change 
of void fraction must be greater than 
some minimum value, in addition to 
satisfying the necessary criterion 
given by equation (9). Further ex- 
perimental evidence is necessary to 
prove or disprove this argument, 
which is based on limited experimen- 
tal data and is presented here more 
as a possibility than a fact. 

Effect of System Pressure 

Similarly, the minimum heat flux 
necessary to cause flow oscillations 
was computed for constant inlet bulk 
temperature by varying the system 
pressure. The results are plotted in 
figure 14 together with experimental 
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Exit saturation temperature, °F 


Figure 14. - Comparison of experimentally determined and theoretically 
predicted necessary heat flux for onset of oscillations as function of 
exit pressure (saturation temperature). Inlet plenum bulk tempera- 
ture, -70° F. 


results for four different velocities. 
Saturation temperature instead of 
pressure is used in the graph for 
convenience. The experimental and 
theoretical results show good corre- 
spondence. It should be noted (see 
table I, p. 8) that, at higher system 
pressure, the oscillations could not 
be obtained in the subcooled regime. 
It may be argued that the increase in 
heat flux necessary to obtain net 
quality at higher pressure is smaller 
than the corresponding decrease in 
the liquid-vapor density ratio. Thus, 
the net effect of pressure increase 
would be expected to be a lower axial 
rate of change of void fraction, which 
is an unfavorable condition for onset 
of flow oscillations in subcooled 
boiling. 

The general trends of the effect 
of system pressure and bulk tem- 
perature found experimentally can 
also be predicted by examining the 
necessary criterion for inception of 
oscillations as given by equation (9). 


With other conditions remaining constant, an increase in system pressure will decrease 


the value of the integral in equation (9) (since t increases with pressure) and thus re- 

s 

suit in a more stable system; a decrease in subcooling will result in a lower temperature 


profile t, which decrease the value of the integral, and thus render a more stable system. 


Effect of Boiler Length 

A comparison between experimentally obtained and predicted heat flux at the onset of 
oscillations with a variable length test section is shown in figure 15. The results are in 
very good agreement and show an inverse proportionality between heat flux and length. 

It should be noted that as the length of the boiler tube decreases the total heat input nec- 
essary for inception of oscillations also decreases, which means that the exit tempera- 
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20 30 40 50 60 70 80 90 100 110 

Ratio of tube length to diameter, L/D, dimensionless 
Figure 15. - Comparison of experimentally determined and theoretically predicted necessary heat flux 
for onset of oscillations as function of boiler length. Bulk velocity, 2. 20 feet per second; inlet 
plenum bulk temperature, 70° F; exit plenum pressure, 20. Opsia. 


ture decreases. In other words, a decrease in length will induce flow oscillations at 
somewhat lower exit bulk temperatures. The results shown in figure 15 were obtained 
at constant inlet bulk temperature and velocity and constant exit pressure. 


CONCLUDING REMARKS 

Since the primary objective of this study was to determine the conditions at the onset 
of flow oscillations, little attention was paid to the amplitude and period of oscillations. 

In general, the amplitude was sufficiently high to allow the instability to be recognized. 
An examination of the traces revealed the frequency of oscillations to be directly propor- 
tional to the system pressure. The velocity, heat flux, and bulk temperature, however, 
seem to have very little, if any, effect on the frequency over the range of this investiga- 
tion. 

A few runs were conducted in which the heat flux was increased beyond the value 
corresponding to the onset of oscillations. It was found that by increasing the heat flux, 
a neutrally oscillating subcooled boiling system could become stable as the thermody- 
namic quality at the exit reached a positive value; increasing the heat flux further would 
cause the system to become unstable again in the net quality region. This phenomenon, 
obtained with an electrically heated tube, is illustrated in figure 16. It should be noted 
that the surface temperature oscillations in the two cases are different in nature. With 
net quality, the temperature is similar to an excursion, usually observed and associated 
with burnout and followed by a slow decrease. In the subcooled case, however, the tem- 
perature rises slowly and then suddenly drops at about the minimum flow, a behavior ob- 
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Figure 16. - Typical flow, pressure, and surface temperature 
behavior in transition region between subcooled boiling and 
net vapor generation. 


served at the incipience of nucleate boiling. 
Some preliminary runs, in which an additional 
thermocouple further upstream was monitored 
on the oscillograph, yielded similar results; 
that is, conditions could be reached, such 
that the upstream surface temperature oscil- 
lated in a fashion characteristic of the incip- 
ience of boiling, while the exit thermocouple 
exhibited the burnout characteristics. Due to 
the oscillatory flow and therefore varying 
local quality, it is not surprising to attain 
this intermittent burnout phenomenon. It is 
believed that the periodic temperature excur- 
sions at the exit are the result of rather than 
the cause of flow oscillations. It is quite 
possible that with insufficient instrumentation 
the quickly developing oscillations can go un- 
detected, resulting in what is usually referred 
to as ’’premature burnout" or ’’submaximum 
critical heat flux. ” 

Some additional runs were made with the 
annular test section where net quality was 
generated at the exit under stable conditions 
by sufficiently raising the inlet bulk temper- 
ature. It was observed that whereas the 
annular flow pattern did exit at the exit, the 
transition from bubbly to annular regimes was 
quite gradual. As the heat flux was increased, 
the transition moved upstream, became more 
abrupt, and finally resulted in oscillatory be- 
havior. This can be explained by the fact 
that in the low quality region the pressure 
drop in the boiler increases with quality, 
with the downstream portion of the boiler con- 
tributing most to the increase. Consequently, 
the local pressure (and hence, the saturation 
temperature) in the upstream portion of the 
test section increased with heat flux. This in 
turn would suppress boiling in the upstream 
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section, which necessitates a higher axial rate of change of void fraction in the down- 
stream part of the boiler. 

It is believed that in net quality boiling, just as in the subcooled case, a critical 
value of thermal nonequilibrium between the two phases must be reached locally (not nec- 
essarily at maximum fluid enthalpy) if the oscillations are to occur. The fact that the 
system remained stable for all air injection tests also indicates this possibility. Further 
studies are necessary, however, to substantiate this theory and determine whether the 
thermal nonequilibrium is the driving force of flow oscillations. 


SUMMARY OF RESULTS 

The criterion given by equation (9) predicted the necessary heat flux for incipience of 
oscillations occurring in low-pressure subcooled boiling systems. The necessary cri- 
terion was derived on the basis of local rates of bubble growth and collapse. The good 
agreement between the theoretical and experimental results indicates that the oscillatory 
flow behavior may be a strong function of thermal nonequilibrium between the liquid and 
vapor phases. 

The experimental observations can be summarized as follows: 

1. Basically, two types of oscillations were observed in subcooled boiling with elec- 
trically heated test sections: constant-amplitude oscillations and oscillations with in- 
creasing amplitude that resulted in surface temperature excursion. It is believed that 
the excursion was the direct result of the relatively low heat capacitance of the tube. 

2. With the aid of high-speed movies, flow patterns at the inception of oscillations 
were studied. It was observed that when flow oscillations were present a quick transi- 
tion from a bubbly type to an annular type of flow pattern existed in the test section. The 
interface oscillated up and down the tube with the same frequency as flow and pressure 
oscillations. 

3. Comparable oscillatory behavior was obtained when boiling was simulated by 
steam injection through a porous wall. When two-phase flow was attained by injecting 
air through the porous wall, however, the system remained stable for all cases investi- 
gated. 

4. An oscillatory, subcooled boiling system could be stabilized by increasing heat 
flux and thus generating net quality. The system remained nonoscillatory until some 
new heat flux value was reached. The surface temperature oscillations in subcooled and 
net quality vapor generation were entirely different in shape. 
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5. With other variables held constant, a decrease in system pressure and an in- 
crease in inlet subcooling tend to destabilize the subcooled boiling system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 1, 1965. 
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APPENDIX A 


SYMBOLS 


c p 

specific heat, Btu/(lb mass)(°F) 

y* 

wall distance parameter, 

D 

tube diameter, ft 


V T iPi 

y, dimensionless 

h f 

fg 

latent heat of evaporation, 


Mi 

Btu/lb mass 

a 

thermal diffusivity, sq ft/hr 

K 

constant 

(3 

heat-transfer parameter, 

k 

thermal conductivity, 


q i 

, dimensionless 


Btu/(sec)(ft)(°F) 


ViVv 5 ! 

L 

tube length, ft 



m 

constant 

e 

time, hr 

N Pr 

Prandtl number, dimensionless 

M 

viscosity, lb mass/(sec)(ft) 

N Re 

Reynolds number, dimensionless 

£ 

coefficient of bubble growth, 

n 

constant 


dimensionless 

P 

pressure, psia 

P 

density, lb mass/cu ft 

q i 

heat flux, Btu/(sec)(sq ft) 

T 

shear stress, lb mass/(sec^ 

R 

bubble radius, ft 

Subscripts: 

r 

tube radius, ft 

av 

average 


V T iPi 

tube radius parameter, r, 

b 

bulk 





Mi 

c 

condenser 


dimensionless 

ex 

exit plenum 

t 

temperature, °F 

i 

inside wall 

t* 

temperature parameter, 

(ti - t)c yfTK 

1- — , dimensionless 

in 

inlet plenum 


SL 

liquid 


q i 

ph 

preheat 

u 

velocity, ft/sec 

s 

liquid-vapor interface 

u* 

velocitv parameter. u 

VWi 

V 

vapor 


dimensionless 

Superscript: 

y 

distance from wall, ft 

C) 

first time derivative 
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APPENDIX B 


UNIVERSAL VELOCITY AND TEMPERATURE DISTRIBUTION 

In order to solve the differential equations for shear stress and heat transfer, 
Deissler assumed that eddy diffusivities for momentum and heat are equal. Furthermore, 
it was postulated that the shear-stress and heat-transfer distributions across the tube or 
boundary layers have a negligible effect on velocity and temperature distributions. The 
molecular shear- stress and heat-transfer terms in the equations were neglected in the 
region distant from the wall. When expressions for eddy diffusivities were assumed, 
Deissler arrived at expressions for nondimens ional parameters for velocity and tempera- 
ture distributions for a constant fluid properties case. The analysis is extended (refs. 8 
and 10) to the case with variable viscosity, a property which for liquids varies more 
rapidly with temperature than other properties. For liquids, the variation of viscosity 
with temperature may usually be represented by 



where m is a constant. When t* and /3 are defined as 


t* 


(tj - 


(Bl) 


and 


/3 = 




equation (Bl) becomes 


— = (1 - ;8t*) m (B2) 

Mi 

When equation (B2) in the shear -stress and heat-transfer equations is used together 
with all the assumptions stated previously, Deissler’ s universal velocity and temperature 
profiles for y* < 26 can be written in integral form as 
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(B3) 


u* 



dy* 

r 


(1 - /3t*) m + n^u*y* < 


1 - 


exp 


n ^ u *y* 

(1 - 0t*) m 


and 



where t*, u*, y*, and /3 are dimensionless temperature, velocity, distance, and heat- 
transfer parameters, respectively, and n is a constant. 

For the region at a distance from the heated wall, that is, for 26 < y* < r?, equa- 
tions containing variable viscosity are neglected and one obtains for universal velocity 
and temperature distributions in this region the following expressions: 

u* = u**_ 26 + 1 In tL (B5) 

y K 26 


and 


t* - ty *_26 + u * " u y*=26 


(B6) 


where K is a constant. 

Equations (B3) and (B4) can be solved simultaneously by iteration. Thus, the solu- 
tion of equations (B3) to (B6) will yield the universal velocity and temperature distribu- 
tions in the tube for any set of conditions. 

The temperature parameter resulting from equations (B4) and (B6) is defined as 


t* 


ftj - t)c p T i 

QiV^i 


(B7) 
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which implies that the shear stress at the wall r. must be known before the temperature 
t can be determined. This can be accomplished by using the 1/7 power law to determine 
the initial value of Tp which was used in the initial value of the heat-transfer parameter 
|3 in equations (B3) and (B4). Simultaneous iteration of these equations yielded the uni- 
versal distributions for the assumed values of r^. From equation (B7) and the definition 
of tg an expression for shear stress at the wall is desired: 



%‘S 


n2 



(B8) 


where 



r r i 

I 1 t*u*(r? - y*)dy* 

Jo \ 

f 1 u*(r? - y*)dy* 

Jo 1 


(B9) 


From equations (B8) and (B9) a new value of Tj is calculated, and it is used for calcu- 
lating the new value of /3. The process is repeated until the calculated and assumed 
values of t ^ correspond within 0. 01 percent. The final values of and t* are then 
used to calculate the temperature profile t from equation (B7). 
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